Constrain on possible pairing symmetry in a two-orbital model of FeAs-based superconductors 



00 

o 

O 
(N 



Wen-Long You,'-^ Shi-Jian Gu,' Guang-Shan Tian,'-^ and Hai-Qing Lin' 

'Department of Physics and ITP, The Chinese University of Hong Kong, Hong Kong, China 

^School of Physics, Peking University, Beijing, 100871, China 

(Dated: July 9, 2008) 

In this work, we establish a few exact identities through commutation of intra-orbital and inter-orbital on- 
site pairings with a two-orbital model describing newly discovered FeAs-based superconductors. Applying 
the conclusion drawn from rigorous relation and physical interpretation, we give constraints on the possible 
symmetries of the superconducting pairing of the model. Hence the favorable pairings in newly discovered 
high-temperature oxypnictide superconductors are proposed. 
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Introduction — The newly discovered family of FeAs-based 
ReOi_xFj:FeAs (Re=La, Nd, Ce, Pr, etc.) high-temperature 
superconductors (SC)[1, 2] spark great interests both experi- 
mentally [3-13] and theoretically [14—32, 34]. So far the tran- 
sition temperature Tc has gone up as high as above 50 K [3- 
5]. The new family material provides another platform to 
explore high-Tc superconductivity besides cuprate supercon- 
ductors. Both FeAs-based materials and cuprate supercon- 
ductors are transition-metal compounds on a two-dimensional 
(2D) square lattice, and their parent compounds are magnet- 
ically ordered. However, there are several significant differ- 
ences in electronic properties. Firstly, the undoped oxypnic- 
tides are bad metals [4, 14], while Cuprates are Mott insu- 
lators. Secondly, neutron scattering experiments have shown 
that the magnetic structure in undoped oxypnictides is not the 
same as a simple antiferromagnetic order in cuprates but in- 
stead a coUinear spin-density wave along the (n, 0) direction 
[4, 9]. Thirdly, probably the most important, multi-orbital na- 
ture of the oxypnictides has been emphasized, in contrast to 
the single-band cuprates. From the band structure point of 
view, it seems likely that all 3d orbitals of the Fe atoms are 
involved in the low energy electronic properties [14, 15]. 

Meanwhile, the pairing symmetry remains controversial. 
Based on fermionic nature of the gap function, the possible 
superconducting order parameters can be classified according 
to the group theory [16-19]. There are surveys that support 
either s- or extended i-wave[20-24],or /7-wave [25, 26], or d- 
wave [27-32],even mixture of Svy and d^i^^i [33, 34]. Never- 
theless, none of them is confirmative. Recently a nodal super- 
conductivity in the electron-doped oxypnictides with multiple 
gaps structure was suggested by specific heat [6] as well as the 
nuclear-magnetic -resonance(NMR) experiments [10-12] and 
point-contact spectroscopy [13]. Therefore, the investigation 
of possible coexistence of various superconducting orders is 
highly desired. 

The Hamiltonian — The Fe atoms in a Fe-As plane form a 
2D square lattice. Due to the buckling of the As atoms, the 
real unit cell contains two Fe atoms. As shown from crys- 
tal field splitting and simple valence counting as well as more 
reliable local-density-approximation(LDA) calculations [14], 
it is reasonable to assume that Fe "idxt and 2idy^ orbitals play 
an important role in the low energy physics of this material. 
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FIG. 1: (color online) Top: The Fermi surface of the two-orbital 
model on the large IFe/cell BZ at yu=1.00, yu=1.45, /j=2.00 from left 
to right. The dashed square indicates the BZ of 2Fe/cell. The param- 
eters are consistent with those in Ref.[35]. 



which are locally degenerate due to the tetragonal symmetry. 
LDA results also show the presence of small Fermi surfaces 
(FS) [23]. In the unfolded Brillouin zone consisting of one 
Fe per unit cell, electron and hole pockets exist around M 
and r points, respectively. The system has multi-band with 
a hole-like Fermi surface around the T point and an electron- 
like around the M point of the Brillouin zone. Upon doping, 
the hole-like FS shrink rapidly, while the electron FS expand 
their areas, as shown in Fig. 1. The dominant contribution is 
clearly around F, M. Thereby, we consider a 2D square lat- 
tice with "dxz, dy" orbitals per site as a starting model to de- 
scribe FeAs-based superconductors [35] for the oxypnictide 
compounds. 
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Here c (d) labels dxz (dyz) orbital, c. (d. ) is creation op- 
erator for electrons of spin o- and orbital d^-r (dy-) at site 
', Hccr = c!^c,>, ni^^cr = dl^dia-, Ti is the pauU matri- 
ces for two orbitals at site i, the hoping integrals are f^. - 

f 1 6{i -i + x) + t26(i -j+y) + h5(i -j + x+y) + ti5(i-j-x'+y). 



TABLE I: One dimensional irreducible representations of Dd, group 
in spatial and orbital space 
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FIG. 2: (color online) The pair coordinates used to defining A_j, 
ranges over 5 positions including (0,0). 



tffa- = t26{i-i+x)+h5{i-i+y)+h6(i-i+x+y)+h6{i-i-x+y) 
^"^^ ^lia ^ ^45(; - i + x-y) + 146(1 - j-x + y)- h6{i - j + 
X + y) - USii - i - X - y) (as illustrated in Fig. 1 of Ref. 
[35]), Uc, Ud are intraband Coulomb repulsion with the rela- 
tion Uc - Ud = U, Ucd interband Coulomb repulsion, and 
Jii the Hund's rule coupling. For later convenience, we set 
fj. - ih + f2)/2, f|| = {h — h)l'^- The space group of Fe atoms 
is P4/nmm, and it is characterized by point group D4 and lat- 
tice translation group T. Therefore, the basis matrix func- 
tions belong to different irreducible representations of point 
group D4, which has five irreducible representations, includ- 
ing four one-dimensional representations (Ai, A2, Bi and B2) 
and one two-dimensional representation (E). Typical bases in 
each representation are listed in Table I. The orbital part of the 
pairing matrix Q. is spanned in the vector space of {djc,, dyS), 
which is an irreducible representation E^ of the point group 
D4, and To,i,2,3 in Table I are transformed as A\, B2, A2, B\ 
respectively [16, 17, 19]. 

In the followings, we will establish sum rules for vari- 
ous pairings by exploiting the commutation relations between 
Hamiltonian and on site pairing operators. Suppose A,, B,, 
C,, and Di be some localized operators defined on lattice A, 
and *?() be an absolute ground state of //a. If they satisfy the 
following commutation relation [//a. A,] - aBj + fiCi + yDj, 
where a ^ 0, jS ^ 0,and y ^ are some constants. Since 
{'i'o\[H^,AiWo} = (Eo-Eo){'i'o\Ai\^o} = so ^(^IB,!^,,) + 
y8<*I'o|C,|*I'o) + r<*I'olAI*I'o> = 0, hence the orders of B„ C„ and 
Di should be either absent together or at least two of them ex- 
ist simultaneously in the ground state [36]. The rigorous proof 
will be stated later 

Intra-orbital pairing — We consider the intra-orbital pair- 
ing. Define on site pairing AJ:"^ = c,.,-fCr,|, ls.f - dr,id,-,i, 
nearest-neighbor (NN) or next nearest-neighbor (NNN) spin 



singlet intra-orbital pairing A*^ 



dr,\d, 



_g- '^'■.T'^,-+a,i~'^''-l'^/-+5,T' 



\dd 



r+S,i 



dr,id t^,and NN or NNN spin singlet inter-orbital 



^li 



pairing A;;^. = c,,td,.^^_^ - c,,id,.^^_^, A^^^. = d,.,^c, 
'^'.i'-/+t?T y calculating the commutation relation with the 
Hamiltonian, which can be expressed as linear combination 
of s and d wave symmetry involving NN and NNN electrons, 
such as on site s wave pairing operator A, = Tjr^r/N, ex- 
tended s* wave pairing operator A,* - 2/(^r-.v + ^r+.t + 
Ar-y + Ar+y)/N, s„, wave pairing operator A.,^^, - 2r(^/-+i+y + 
Ar-A-v + Ar+.t-v + Ar-_uy)/N, d^^-yi wavc pairing operator 
Aj,2_,2 = i;r(Ar-.; + Ar+.;-A,._v-A,.+v)/A^, and £/.j7 wave pairing 
operator A^^,,, = Y^A^r+m- + K-x-y - K+x-i - K-x+y^N. The 



respective pair coordinates are shown in Fig. 2. 
Thus, we have. 



[H,AT] 


= t^AZ + hA^_^ + t.AZ, 
- (U-2fi + 2Ucdnr,d)A7 , 


- '*< 


[H,Af] 


* \dd . \dd , ^ \dd 

= t±A„ - t\\A^^^_^^ + hA,^^ 
- (U-2fi + 2Ucdnr,c)Af . 


-t<. 



(4) 



(5) 



Since there are inter-orbital hopping in the original Hamilto- 
nian (1), inter-orbital pairings appear on the right hand side 
of (4) and (5). Furthermore, for quartic terms such as iirjAf, 
we treat them in the mean field sense by taking number oper- 
ators n,. as its average (n,), assuming that there is no charge- 
density wave (CDW). Even though CDW occurs, the strong 
Coulomb repulsion will destroy on site s wave pairing, i.e., 
A5=0, and it will not affect our conclusion. Thus, in the fol- 
lowings, we treat the number operators for two orbitals as 
(«r,d) = (ir,e) = («) in all quartic terms. Make a linear combi- 
nation of (4) and (5), 

[h, Af - Af ] = tAAZ - Ait) + f||(A- _^^ + Aif^J 
+f3(A^:;,. - Af ) - (U-2fi + 2Ued{n})(A:'-A'/). (6) 

To express the form more explicitly, with the help of 
Nambu representation (pa-{k) = ic-k,a-,d-k,cr), the pairing 
gap in momentum and orbital space can be expressed as 
I^kafi 'l>oA-k)f(k)iTi)a/}(f>/},cTil^)' whcrc f{k) is symmetry fac- 
tor in momentum space, 

[H, T3] ~ fx(cos kx + cos ky)T3 + f||(cos kjc - cos ky)To 

+ ty,{cos kx COS ky)T3 - (U -2fj. + 2Ucd{n})T3 (7) 

There are four types of pairing patterns on the right side of 
equalities (7), listed as No. 1-4 in Table II. 

Inter-orbital pairing — When the inter-orbital pairing is 
taken into consideration, more pairings are involved. Define 
on site inter-orbital pairing as AJi"' - Cr^idr^i, A'f' = dr,^Cr,i, 
and NN or NNN inter-orbital pairing operators in a similar 
manner, 

[h, Af + Af ] = f^(A^f + Af ) + f3(A^f^, + Af ) - t4{A''dl + ^t) 
3Jh\, , 3Jh 



[2p-U- 2Ucd + -^){n) ■ 



(Af + Af ) 



(8) 



The pairing gap in momentum and orbital space, 

[//, Ti] ~ t4^(CQSkx + CQSky)T\ 



TABLE II: Potential pairing basis matrices under different irre- 
ducible representations of the model, which are classified into three 
groups in the table. The first column is the index number, the sec- 
ond and the third columns denotes the representations and the ba- 
sis matrix functions respectively. The parities of spins [singlet(S) or 
triplet(T)] and orbitals [symmetric(s) or antisymmetric(a)] are shown 
in the forth and the last columns respectively. 



No. 


IR 


Basis 


Spin 


Orbital 


1 


fii 


(cos fcj- -1- cos A;,.)t3 


S 


s 


2 


fii 


(cos ky: - COS A;,,)to 


S 


s 


3 


fii 


(cos ky: COS A;,,)T3 


S 


s 


4 


fii 


T3 


S 
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fi2 


(cos kx + COS ky)Ti 


S 
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fi2 


(cos k^ cos ky)Ti 


S 
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Bi 


{sinky smky)To 


S 
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Bi 


Tl 


S 


s 
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Ai 


(cos ky + cos ky)iT2 


T 


a 


10 


Ai 


(cos ky COS ky)iT2 


T 


a 


11 


A2 


iT2 


T 


a 



-I- f3(cos kj, COS ky)Ti - f4(sin kj, sin ky)To 



(2^-f/-2t/., + H^)<«>-H^ 



Similarly, with the definition of spin triplet pairings Ai^ 

6 

'^'■T^.+ai + crid^^g^, Af = drjc^^g^ + t/^c^^g^, we have 

[//,Af-Af] = fi(Ajf - At") + f3(A^'f, - A^;;,) 

Jh\, ., Jh 
— ]{n) + — 



[2^i-U - 2Ucd-^Yn} + ^ 



(Af-Af)(10) 



[H,iT2] ~ t j_(cos kjc + cos ky)iT2 + t^icoskjc cos k^)iT2 

IT2 (11) 



{2^-U- lUcd - ^-f){n) + ^ 

There are four types of pairing on the right side of equahties 
(9) Usted as No. 5-8 and three in (11) Usted as No. 9-11 in 
Table II, where the gap behaviors are also enumerated. 

Inequality and analysis — With the use of the above com- 
mutation relations, let us establish the sufficient condition 
for the coexistence of two long-range orders rigorously. We 
should take advantage of off-diagonal long range theory [37] 
and generalizing the approach of Ref. [38] to obtain a strict 
proof. 

For a general operator Gk defined on lattice A, it has a long- 
range order if and only if (iff) its reduced-density matrix in k 
space satisfies 

<^o(A)|G|G,|1'o(A)> >ANty, A>Q, 

at a certain k - ko point as A^a — > oo with fixed density. As- 
suming an operator At satisfies the following commutation re- 
lation 



[Ha, At] = aBk +/3Ck + yD^ = Qk , 



then we can prove that 

<^o(A)ie^ei|To(A)) < m(Qk)m(Ak) , (13) 

m(Qk) = 7<**'o(A)|[e^, [Ha,Q,]Wq(A)} , (14) 

m(At) = ^<^o(A)|[A^,[//a,AJ]|>I'o(A)) . (15) 

Here m(Qk) and m{Ak) are quantities of order (9(1) as A^a tends 
to infinity. Therefore, the correlation function of Qi is at most 
a quantity of order (9(1). Then expanding the inequality (13) 
and applying the Cauchy-Schwarz inequality, we have 

^ |^^|2<>I'o(A)|Gf G^|^o(A)> < m(Qu)m{A,) 
+ ^ 2\A^Uy\ ^<To(A)|Gf' G^|>I'o(A))<*Po(A)|G;"'"G^|>I'o(A)) 

where G^ represents Bk{Ck,Dk) in Eq. (12). Now let us as- 
sume that *l'o(A) has a kind of long-range order, say Bk, then 



<>Po(A)|B ,Bi|^„(A)) > ^bA^a, ifi > , 



(16) 



as A^A — » °°- Since the left hand side of inequality (16) is 
n . (9) a quantity of order 0(Na), and other operator's correlation 
function is, at most, a quantity of order (9(1) in the thermo- 
dynamic limit, the right hand side of the inequality can be, 
at most, a quantity of order (9( VAa). Therefore, *I'o(A) must 
have another long-range order of, say Ck, that is, 

(%(A)iq^C,|>Po(A)) > AcNa, Ac > 0. (17) 

The conclusion is that either all orders are absent or at least 
two long-range orders must be present simultaneously in the 
ground state *1'()(A) of the Hamiltonian Ha- 

Applying the above rigorous result to Eq. (1), we have a 
basic conclusion that orders of the same group listed in Table 
II, say No. 1-4, should either are all absent, or at least two of 
them coexist. Specifically, we can refine this conclusion as 
follows. 

Firstly, the strong on-site Coulomb interaction and large 
Hund's coupling suppresses all the states involving s wave 
pairing, then No. 4 (symmetric intra-orbital s), 8 (symmet- 
ric inter-orbital s), 1 1 (anti-symmetric inter-orbital s) will be 
unfavorable. 

Secondly, Eq. (7) shows that all order parameters in the 
first group of Table II are of spin singlet, intra-orbital, and 
even parity pairing symmetry. Nevertheless, our rigorous re- 
sults [Eqs. (16 and 17)] require that orders of No 1-3 should 
either vanish simultaneously or at least two of them coexist. If 
two of them coexist, the magnitudes of the two order parame- 
ters are determined by the multiplying their pairing functions 
(/i (k), f2{k)). So if s* and s^y coexist, /i (k) = (cos k^ + cos ky), 
f-i{k) - cos kx cos ky. Their overlap becomes dominated 
around the hole pocket about F points in the Brillouin zone, 
but is very small around the electron pocket about M points. 
(12) While if d^i-yi and Sxy coexist, f2{k) - (cosk^ - cosA:^,)^ 



fsik) = cos ^t cos fey, their product has an enhanced contri- 
bution from the electron pocket about M points, but is sup- 
pressed from hole pocket about F points. However, the coexis- 
tence of s* and dj^i^yi is not favorable because of their overlap 
is very tiny around both M and T points. 

Thirdly, in the second group of Table II, there are No. 5 
(symmetric inter-orbital s*), 6 (symmetric inter-orbital s^y), 
and 7 (intra-orbital dxv) left on the right side of (9). Since 
fi(k) = sin ^x sin ^y is peaked around (nl2,nl2) , if No. 7 co- 
exists with No. 5 or No. 6, fi{k) has a tiny contribution at M 
and r points and does support the current fermi surface topol- 
ogy, therefore No. 7 is not favored. So the remaining s* and 
Sn- of spin singlet must be either absent together or coexist 
simultaneously. According to the behavior of symmetry fac- 
tors fs((,)(k), the coexistence of s* and Sxy is competitive only 
around hole pocket about F points, but not favored in electrons 
pocket about M points. 

Fourthly, regarding with the third group including No. 9 
(antisymmetric inter-orbital s*) and No. 10 (antisymmetric 
inter-orbital Sxy) on the right hand side of (11), our rigorous 
results [Eqs. (16) and (17)] impose that either both of them 
are absent or coexist. Because both of them carry antisym- 
metric orbital parity and become gapless in excitation spec- 
trum, which are inconsistent with the experimental evidence 
of nodal gap [12, 13], the chance of coexistence seems slim. 

Finally, in the weak coupling limit, two orbitals' energy 
splitting might lead to a mismatch of inter-orbital pairing in 
momentum space with opposite sign, instead of pairing be- 
tween two different |A:|s. That is the piling up of low-energy 
density of state in the gapless SC state will lead to a Fulde- 
Ferrel-Larkin-Ovchinnikov state with magnetic ordering, and 
does not create SC instability. In this sense, the orbital anti- 
symmetric pairing state such as Nos. 9 and 10 might be ruled 
out. Moreover, other inter-band pairings, such as Nos. 5 and 
6 in Table II are also not favorable according to the analysis 
based on the FFLO state. Then we arrive at our further con- 
clusion that around half filling, in the electron doping region, 
the system will favor coexistence of dj2-y2 and Sxy waves pair- 
ing, while in the hole doping region, the system might prefer 
to have s* and Sxy waves pairing. 

To summarize, we have built some identities based on a 
two-orbital model, and obtained constraints on a few possi- 
ble pairings. Our results provide more information than group 
theory classification. According to the sufficient condition for 
coexistence of two superconducting orders and resorting to 
physical consideration, we propose the most favorable pair- 
ings around half filling. Although our discussions are based 
on a two-orbital model, it is straightforward to generalize the 
strategy to other Hamiltonians even if more orbitals are in- 
volved. In principle, we have not ruled out the spatial odd 
parity pairing, e.g., p wave, which can be achieved by com- 
mutation between an odd parity pairing operator and Hamil- 
tonian, and the sufficient condition of coexistence of the odd 
parity pairings is still applicable. Nevertheless, they do not 
get along well with our fermi surface topology analysis given 
above. 
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